NASA-CR-197671 


/y/Uy 

030601-4-T 


Simulation of Thin Slot Spirals and Dual 
Circular Patch Antennas Using the Finite 
Element Method With Mixed Elements 


J. Gong, J.L. Volakis and M.W. Nurnberger 


i 


National Aeronautics and 
Space Administration 
Langley Research Center 
Hampton, VA 23681-0001 



January 1995 

THE UNIVERSITY OF MICHIGAN 

Radiation Laboratory 
Department of Electrical Engineering 
and Computer Science 
Ann Arbor, Michigan 48109-2122 
USA 


0 I o 

f*- | oo >t 

O" 3 O o 

ro rn 

1 X I — o 

iThir g rn 
o- \ o c o 

7 i Z D O 


rsl 

rn 



< 00 v ++ 

-J H- O i 

3 2 * > 

UL Vj W UJ W 

O OC H X CD C 

h m uj 

z u z j 
O UJ * c 

m jii 4-> m 

< Q W ^ 

< D UJ UJ O ■*- 

-jo xx ax 

3 K i-i a? u 

X O X C* - 

*-» Z o X 

00 < z: x ^ 

oo un ««* a? 

^ j 3 3 l ir 

<t cn o 
a£ oo O 0 0" 

^ h < o u ^ 

Nazia 

0 s tA Z b* • 

rH UJ UJ C 

1 H K X 03 <*0 

a: o -z 3 n 

u j < c 

I 00 2 T c cn 

< X tu fC C 

z g i - a 

< H h III E 3 

Z X < -J 0> C rJ 

w v- a uu oo 0) 





030601-4-T 


NASA Grant No. NAG-1-1478 


Grant Title : 


Simulation of Conformal Spiral Slot Antennas on Composite 
Platforms 


Report Authors : J. Gong, J.L. Volakis and M.W. Numberger 

Report Title: Simulation of Thin Slot Spirals and Dual Circular Patch Antennas 

Using the Finite Element Method With Mixed Elements 


University Project Director : John L. Volakis 

volakis@umich.edu 
Telephone: (313)746-0500 


NASA-Langlev Project 

Monitor : Fred Beck 

Telephone: (804) 864-1829 


University Address : Radiation Laboratory 

Department of Electrical Engineering 
and Computer Science 
The University of Michigan 
Ann Arbor, MI 48109-2122 


Date : 


January, 1995 


l 




Table of Contents 


Page 


Preface 

1 . Slot Spiral Antenna Modeling Using Hybrid/Mixed 

FE-BI Technique 

2. Design Considerations for an Archimedean Slot 

Spiral Antenna 

3 . An Efficient and Accurate Model of the Coax Cable 

Feeding Structure for FEM Simulations 

4 . Arbitrarily Shaped Dual-Stacked Patch Antennas: 

A Hybrid FEM Simulation 

5 . Performance of an Artificial Absorber for Truncating 

FEM Meshes 


v 


1 ^ 
■ • A 

21 o 
47) 



77 0 


m 


P PAGE BLANK NOT FILMED 

I* 

PAGE Li iNTEft ncftou . . .m 




IV 



C/W/S 


Preface 


This semi-annual report describes our progress up to mid-January 1995. The report contains five 
sections all dealing with the modeling of spiral and patch antennas recessed in metallic platforms. 
Of significance is the development of decomposition schemes which separate the different regions 
of the antenna volume. These schemes are used in modeling the spiral and the dual circular patch 
antennas. In the case of the spiral, a mixed (edge and node) element formulation is developed to 
avoid use of small elements in modeling the slot. This also allows for use of the simpler wire 
green's function in modeling the radiated fields. Overall, the resulting modeling approach leads to 
substantial simplifications and is perhaps the most efficient within the context of hybrid FEM 
formulations. The implementation of this formulation is nearly complete but we already fabricated 
and tested an archimedean slot spiral operating at a center frequency of 1.1 GHz (see Section 2 of 
the report). The main goal of the experiment was to examine the feasibility of a new feed design 
for the slot spiral antenna. This feed design consisted of a microstrip line below the metallic 
section of the slot spiral but, unfortunately, our first design trial required improvements to bring 
the axial ratio close to unity. A new design is currently under fabrication and with the help of the 
new software, it is our hope that the feed design issue will be resolved over the next two months. 

Substantial effort was devoted over the past six months in improving the feed model in the 
context of the FEM. This is addressed in Sections 3 and 4 where we describe simple, but very 
successful implementations of the coax probe and slot-coupled microstrip line feed. In the latter 
case, the feed network is modeled using brick elements whereas the antenna volume is modeled 
using tetrahedrals. The two computational regions were then coupled by enforcing potential 
continuity (a new approach) across the slot joining the two regions. A design of a dual stacked 
circular patch configuration is included which demonstrates the validity of the potential continuity 
condition. Also, as noted in Section 3, the proposed coax cable model yields input impedance 
results which overlay the measured data. 

Finally, Section 5 of the report presents an innovative scheme for truncating finite element 
meshes. This is an alternative to the absorbing boundary conditions and employs a perfectly 
matched anisotropic absorber. It was already tested for terminating transmission lines in 
microwave circuits and we are in the process of implementing it in our FEM antenna codes. 

In the next few months our goal is to complete the spiral antenna code and employ it to examine 
several spiral antenna configurations operating at different bands, including that from 944- 
1256 MHz as discussed in the proposal. The code will be also enhanced with several features 
relating to feed modeling and mesh truncation. Finally, additional fabrications and antenna 
measurements will be made. 
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Slot Spiral Antenna Modeling Using 
Hybrid/Mixed FE-BI Technique 

Jian Gong and John L. Volakis 
Radiation Laboratory 

Department of Electrical Engineering & Computer Science 

University of Michigan 
Ann Arbor, MI 48109-2122 


Introduction 

This part of the report is concerned with the numerical simulation of the , 
printed Slot Spiral Antenna (SSA) using the finite element - boundary in- 
tegral (FE-BI) analysis. It has been reported that the FE-BI technique is 
suited for modeling microstrip patch antennas of any shape, printed on lay- 
ered planar structure or cylindrical platform and fed with a coaxial cable or a 
microstripline network underneath the radiating elements (the network feed 
modeling is referred to the attached article). However, direct application of 
the hybrid FE-BI technique to thin slot spirals requires excessive sampling 
rates to accurately simulate the geometry. This is particularly so when edge 
elements are used for tessellating the cavity volume. Although the edge ele- 
ments are better suited for modeling metallic edges and corners they cause 
major meshing difficulties in the case of thin slots in a metallic plane. The 
edge elements must be specifically designed to have one of their edges per- 
pendicular to the slot and this places severe restriction on the volume mesh in 
addition to increasing the computational requirements. Moreover, the FEM 
computational domain beneath the slot spiral antenna may contain a feed- 
ing circuitry or dielectric substrate whose structure may have a significant 
difference in size than that of the thin slot spiral and as a result a large, 
ill-conditioned system may be obtained. 

To alleviate the meshing/modeling difficulties encountered with cavity- 
backed slot antennas, in this report we describe a mixed finite element- 
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boundary integral formulation. As in the past, the boundary integral is used 
to describe the radiation of the slot above the cavity. However the cavity is 
now modeled using a suitable mix of edge and node elements. The latter are 
used only at the aperture of the thin slot so that the nodes follow the center 
line of the slot. In this manner, regular size elements can be used regardless 
of the slot’s width and any meshing restrictions are substantially relaxed. 
The proposed mixed element FE-BI formulation introduces three different 
computational regions and as expected this complicates the generation of the 
discrete system. In the following sections, we describe the proposed mixed 
element FE-BI formulation and proceed with the development of the system 
by ensuring field continuity at the interfaces of the computational regions. 
We conclude with the presentation of some preliminary results. 


Mixed FEM Formulation 


Consider the configuration where a thin slot spiral antenna (SSA) is etched 
on a dielectric substrate recessed in a cavity enclosed by perfectly electric 
conducting (PEC) walls. For the purpose of this analysis, the cavity region 
will be considered as inhomogeneous and of arbitrary geometry. To solve for 
the field in the interior and exterior region via the finite element method, 
standard approach is to minimize the functional 

JF= - f {(V x E)— • (V x E) - kltrE • e} dV + jk 0 Z 0 [ E x H -ndS, (1) 
2 JV l ) JS 

where S denotes the slots cross-section, V is the cavity volume below the slot 
and (E, H) are the electric and magnetic fields, respectively. In discretizing T 
we must introduce appropriate field expressions in each of the computational 
regions. On the slots surface the boundary integral will be used to relate the 
tangential components of E and H. This yields a functional only in terms of 
E, and upon differentiating T with respect to E we then obtain a system of 
equations for the solution of the node/edge fields. Since the cavity volume 
consists of node and edge elements, we consider two field expressions. For 
the eth element these are of the form 


/ ELi E k V k ( r) r € Vi 

\ Ej=i (*E xq + yE yq + zE zq )N q { r) r € V 2 


( 2 ) 


where Vi and V 2 denote the volumes occupied by the edge and node ele- 
ments, respectively. The coefficients Ek represents the edge fields at the kth 
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edge-based shape functions defined in [1]. Similarly, A, denotes the node 
based shape function for the node q and E xq is the field coefficient for the 
x-component of the nodal field. The node— based shape functions are given 

by 


where 
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(4) 


(5) 


in which the indices i,j,m,n are in the cyclic order and represent the node 
numbers of the eth element. 

A special attention shall be paid to the treatment for the nodes in the slot, 
along whose axis the tangential E component must be set to zero in order to 
satisfy the boundary condition. This can be accomplished by rewriting the 
second expression in (2) as 


E e = ^ (xE xq + i/E yq + zE zq )N q (r) + {rhE m4 + zE Z 4 )N 4 (t ), (6) 

9= i 

if the 4 th node happens to be in the slot, and rh — z x l. It is observed that 
this boundary condition treatment indeed reduces the degree of freedom for 
each node from three to two. The nodes on pec walls may similarly be dealt 
with using this technique and in particular the degrees of freedom for each 
node on pec walls becomes one. Therefore, one has the expansion 

E e = Y, ( xE xq + yE yq + zE zq )N q {r) + zE z4 N 4 ( r), (7) 

9=1 
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assuming the 4 th node is on a pec wall. If an element's facet coincides with 
the pec ground plane, the degrees of freedom associated with this element will 
be 6 rather than 12, as is usually the case. The implementation of boundary 
condition is treated at the stage of entry generation and matrix assembly to 
avoid redundant matrix size declaration. 

From the above, the field everywhere in V = V x + V 2 can then be expressed 
in discrete form as 


E = EE £ , 
e 


(8) 


and when this is substituted in (1) we obtain the system of equations 


£ K’ 


{Ek] + X [a { 2) ] {E qx } + {B I - subsystem] = 0. 

C2 6 V 2 


(9) 


These element equations must be assembled to generate a partly sparse, 
partly full system of equation, and in doing so it is necessary to identify 
the relationships among the edges, nodes and boundary field coefficients. In 
addition, the node shape functions defined in (3) are not divergence free, 
i.e. they do not satisfy the V • E = 0 condition, and thus to eliminate 
non-Maxwellian solutions, it is instructive to add the term 



(10) 


to the matrix element of [A (2) ] . Note the sum in {} is over the u (= x,y and 
z) components. The assembly of the element equations is more challenging 
for our specific case because of the mixed volume elements and their interface 
with the boundary elements. Before addressing this issue, we first describe 
the formulation of the BI subsystem, thus completing the specification of the 
element equations (9). 


Boundary Integral For Thin Slot/Wire Struc- 
ture 

As already presented in the above section, the functional (1) contains a term 
involving the tangential component of the magnetic field, h x H, on the 
aperture S or the slot area(s), where n is normal to the aperture and for 
our case n = z. This field should be represented via an integral equation in 
terms of the electric field E or the equivalent magnetic current in the thin 
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slot. To construct the integral equation, we start, with a mixed potential 
integral equation for thin wires of arbitrary shape. Figure 2 shows the top 
view of a thin spiral slot antenna which is located on the boundary separating 
the exterior free space from the finite element 3D domain (cavity). 

The complementary physical problem corresponds to an electrically con- 
ducting thin wire of arbitrary shape carrying the electric current I e along 
only the axial direction. The associated integral equation can be written as 


/ • {E = — juA — V</> | 


/ • E‘ r <E source gap 
0 otherwise 


(11) 


where E‘ represents a plane wave incidence or a delta gap source, l is a unit 
vector tangential to the thin wire at point r and 


A = g 



dl, 


, .V-A 

<P = 3 

up 0 eo 


-^-V • [ Gle dl. 
U>fi 0 Jl 


(12) 


In these, I e is the electric current and G the free space Green’s function 


G{ r, r') 


e -i*|r-r'| 

AttloqCo 


On substituting (12) into (11), we obtain the integral euqtion 


= i-jkz 0 J i [ 


*+F vv l 


G ■ Ie dl + 


i • e 1 ' 

0 


r € source gap 
otherwise 


( 13 ) 


in which 8 is the unit dyad. 

The quadratic basis functions 

r(0 = 2(£ - -)(^ - l)r n _i + 4^(1 - £)r n _x/ 2 + 2^(^ — 2) Fn ’ (^) 

where n = 1,2 , . . . , iV, are chosen to expand the electric current I e and £ 
is the simplex coordinate whose range is [0, 1]. Galerkin technique for the 
residual testing is employed for the method of moment (MoM) implementa- 
tion of (13), and it should be remarked that the thin structure assumption 
(which is adequate for spiral antenna analysis) leads to a one dimensional 
boundary integral subsystem. This significantly increases the efficiency of the 
simulation and the geometry is prescribed using an independent pre-process 
program to generate a ffive— point per record’ (corresponding to two adjacent 
segments) data set as an input file. The singularity for the self-cell entries 
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Figure 1 : A typical spiral slot etched on top of ground plane. The spiral is 
defined as r = a& + b for one arm and a, b may be determined by a = 

6 = T\ — a6 1 = r 2 — a$2, where r lj 2 0i,2 represent the polar coordinates for the 
inner ends and outer ends, respectively. 




is dealt with by assuming a finite width spiral slot structure and this leads 
to a simple numerical evaluation of the matrix elements. The accuracy of 
this implementation was validated by comparing with the referenced data 
available in the literature [2, 3, 4]. Figure 2 and 3 show the comparison of 
real and imaginary parts for the electric currents generated by a voltage gap 
excitation on a circular wire ring. Indded an excellent agreement is observed. 

Additional results for different thin spiral configurations are given later in 
this part of the report and these clearly show the promising capability of the 
technique. 

Returning to the surface integral of then functional (1), we have 

jkZ 0 J s E x H-hdS = -jkZ 0 j(E x z)W -Udl (15) 

with dS = W dl and W is the width of the slot. If W(E x z) is defined as 
an equivalent (line) magnetic current, the magnetic field H measured at the 
z = 0 plane can be expressed as 

H = -jk 0 Y 0 J i G m -I' m dr (16) 

where G m is 

5„ = 2|s + -4vv)g, 

in which the factor of 2 accounts for the ground plane contribution. Then 
(15) becomes 

—jkZ 0 /(E x z)W ■ H dl = -k 2 Q jjl m ■ G m -\ m dl' dl (17) 

Note that following the standard discretization procedure, the magnetic cur- 
rents I m and Y m in (18) are expanded using the quadratic basis functions 
in (14) in terms of a one dimensional variable £ and this variable number 
reduction is due to the application of a prescribed geometry date set, as 
aforemention. We observe that the unknown coefficients for both currents at 
the observation and source points, respectively, should be treated in different 
manner. That is, the differentiation for extremization is only applied to the 
current coefficients at the observation points, remaining those at the source 
points constant. This formulation leads to a same discrete form integral 
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Imaginary Part of Current (Am) 


Imaginary Part of Current on Loop Wire 



00123456789 

Segments 


Figure 3: Imaginary part of current in the same circular ring wire antenna 
as in figure 2. 
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equation as the one derived from (13), where Galerkin's technique was used, 
except for a constant factor. 

It should be noted that in implementing the method of moment for the 
slot spiral, the finite width W of the slot has to be involved (to maintain the 
correct units) and indeed it is an important geometric parameter for the spiral 
description when the Green's function is evaluated at the self- cell terms. 
In addition, the forms for the dyadic Green’s function and the equivalent 
magnetic current are not arbitrarily chosen; usually one determines the other. 
However, the overall result should always give rise to the same answer. 


System Combination 

We describe in this section the issues relating to the interface of edge and 
node-based element regions, as well as the generation of the overall system 
for solution. For discussion purposes, the following notation is used in this 
section: 

• Three subsystems: 

a. edge-basis FEM: E° 

b. node-basis FEM: E 6 

c. BI: E s 

• Basis functions: 

a. edge-basis: V, 

b. node-basis: Nj 

c. BI: S q 

It is known that the hybrid FE-BI system is conventionally combined 
by enforcing the continuity of the tangential electric fields on the aperture 
linking the cavity field unknowns and those in the boundary integral. For 
patch antennas, since the aperture size is relatively large around the patch, 
this practice is appropriate. However, in the presence of the thin slots, this 
scheme is no longer valid because it would lead to very fine meshes and thus 
a large number of unknowns. This difficulty is avoided herewidth by making 
use of nodal elements in the thin slot and its adjacent region. In particu- 
lar, applying (6) and (7) for field expansions offers the inherent opportunity 
of combining the fields above and below the slot without a need to either 
increase the mesh sampling rate or carry out interpolations, both of which 
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worsen the system condition. The way of combining the node based FEM 
system b with a thin slot BI is similar to that for patch antenna simulation 
where the cavity unknowns beneath the ground plane coincide with the BI 
unknowns. Here, we choose nodal basis elements and the field unknowns in 
the slot are defined on specific nodes in both system b and the BI subsystem. 
A one-to-one correspondence of the fields in these two subsystems allows a 
direct combination (see [1]). 

Consider the two separate FEM subsystems (one for nodal elements and 
the other for edge elements). It is apparent that explicit coupling across their 
interface is needed and in this report we propose a technique to couple the 
node and edge-based FEM systems without altering the condition number. 
Consider an interface T between region Vi (edge elements) and V 2 (node 
elements). Since the tangential field continuity requires 

t-E {a) = i -E (t) , rer (18) 

where i is a tangential unit vector in T and 

E (a) = £ 4 0) S,(r), E (6) = ]T EfN^r). (19) 

t 3 

The expansion functions S,(r) lie in the interface T, i.e., there is no normal 
component perpendicular to the interface; but functions Nj(r), in general, 
have both components. If m is defined to be normal to the interface, then 
(18) may be rewritten as 

m x m x E (6) = E (a) r e Y (20) 

or 

m x m x £ Ef N,-(r) - £ JS?{ a) Si(r) 

3 i 

= ££f¥ J (r)-££|->S,(r)=0 (21) 

3 i 

where Pj(r) = m x m x Nj(r). We introduce a vector residue 

R = Y.E? ) ri(r)-'LE^S i (r), (22) 

3 i 

and a functional term for coupling the system a and b may then be given by 

T ah = Jx-RdS. (23) 
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A is called a vector Lagrange multiplier (as opposed to a conventional scalar 
multiplier) and it may be expanded in terms of vector basis functions in 
the same fashion for vector fields E. In particular, we assume an expansion 
function T p and 

A = £ A p T p (r). (24) 


Hence, 


I ^ ApTp(r) ■ (£3 j( r ) — £ E^S,(r))dS. (25) 

,/r P 3 i 

Differentiation with respect to A p , E, and Ej yields, 

/ T p • Pj DS — ^2 ^«' a) / T p • S, DS 

? jr i jr 

£ A p / T p ■ P, DS 

v Jr 

— XI Ap / Tp • Si DS (26) 

p Jr 

These are the explicit expressions for coupling the subsystems a and b (first 
two discrete terms in (9)). It is observed that when the expansion function T p 
is proporly chosen, the entire discrete system of equations is symmetric after 
combination of the three subsystems. This is one of the most significant and 
attractive features of the formulation proposed and described in this report. 

It should be remarked that this hybrid FEM formulation, regardless of 
its complexity, is promising and capable of simulating a thin slot antenna 
supported by a cavity or cascaded cavity structure that is recessed in a planar 
or curved platform. More composite configurations such as aperture coupled 
feeding network or inhomogeneous dielectric fillings may be tackled with no 
extra burden. 


dE ab 

dX p 

d? ab 

dEf 

dT ab 

sip 
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Figure 4: Geometry of a two arm Archemedeam Spiral in refernece [4]. The 
line connecting the outer end to the center is not present in practice and the 
feeds are located at two inner ends. 
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Figure 5: Predicted current distributions on the spiral in figure 4 and an 
excellent agreement is observed in comparison to the results shown in [4] . 
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Figure 6: Axial ratio vs. 9 for the two-arm spiral antenna shown in fihure 1. 
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Figure 9: Radiation pattern for the spiral in figure 1. 
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Figure 10: Admittance vs frequency for the spiral antenna in figure 1. 
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Design Considerations for an Archimedean 
Slot Spiral Antenna 


1 Design Goals 

Develop a 118-157 MHz, vertically polarized, low-profile (or conformal) an- 
tenna as a replacement for VHF AM blade antennas on aircraft. This design 
is to be arrived at by scaling the dimensions of an antenna designed for a 
center frequency of 1.1 GHz. The design prior to scaling may have the follow- 
ing maximum dimensions: Diameter < 3.70" and Thickness < 0.50". The 
suggested/required antenna type is a 4-arm spiral antenna. 


2 Initial Comments 

Although a 4-arm spiral design was originally suggested, a 2-arm spiral may 
also be used, as both mode-1 and mode-2 (sum and difference) radiation 
patterns aren’t required. While a 4-arm spiral can easily be designed should 
both sum and difference patterns be required, the 2-arm design will provide 
the required sum pattern and simplify the design problem somewhat— only 
one feed is required, and the feed area geometry is more straightforward. 
Polarization requirements dictate that a slot spiral be used, as opposed to a 
wire spiral. 


3 Potential Antenna Designs 

Two similar radiating structures were considered. The first is the standard 
archimedean spiral antenna, as shown in Figure la. The second is a hollow 
archimedean spiral antenna — essentially a standard archimedean spiral with 
the inner portion removed, depicted in Figure lb. The design shown in 
Figure la is known to work, while the modification shown in Figure lb offers 
the potential of an integrated feed network and a fully planar design. This 
will be discussed in more depth in a later section. 
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4 Potential Feed Designs 

As mentioned above, two different antenna structures were investigated Ge- 
ometrically, they are quite similar, but they allow two fundamentally different 

fee The 4 S fir n st class of feeds make use of the standard archimedean spiral ge- 
ometry shown in Figure la. While the combination of this antenna and these 
feeds extends the bandwidth of the system beyond that indicated in the g 
inal specification, it also allows the use of standard spiral antenna design and 

feeding techniques. _ ~ 

Several standard feeding techniques exist for spiral antennas. One com- 
mon method uses a balanced two-conductor transmission line connec ed on 
the underside of the antenna, across the slot at the center of the spiral [2 8 . 
Another is the so-called “infinite-balun” coaxial feed, where a coaxial cable 
spirals in toward the center of the spiral along the underside of one of 1 1 
metallic arms, and is connected across the slot at the antenna s . center Fre- 
quently a dummy cable is added along the underside of the other metallic 
arm to keep the structure symmetric and the impedance consisten l , , J- 
The “infinite balun” feed has the advantage over the balanced two- wire 
feed for simple structures - the balun is inherent in the feed design, is as 
broad-band as the antenna, and is simple to build. More complex systems, 
such as multi-mode direction-finding systems, etc., often use some variation 
of the balanced twin-wire feed, where the balun can be integrated moo e 
parts of the system, such as the beamforming circuitry [8]. When accuracy is 
not as important, the balun may also be placed inside the absorbing cavity 

[4] ' The second type of feed, shown in Figure 2, makes use of the unused area 

in the middle of the hollow spiral shown in Figure lb by placing the balun 
structure there. Essentially, a coaxial line is brought up through the cavity, 
as the balanced twin-wire line would be for the first feed and transitions to a 
co-planar waveguide (cpw) transmission line on the top of the substrate. Thi 
cpw line tees, and the two resulting cpw lines are then converted to slotlmes 
through balanced cpw-slotline transitions [5,6], and become the radiating 
arms of the slot spiral. 
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5 Discussion of New Antenna and Feed De- 


signs 

• Coax-cpw-slotline Feed 

As mentioned in the above description, this feed, pictured in Figure 2, 
makes use of two cpw-to-slotline transitions, which also act as baluns. 
These transitions are fed by a coax-to-cpw transition, allowing a single 
coaxial feed. Initially, not much consideration was given to the over- 
all size of each of the baluns. It was expected that, although little 
room would be available, the feed could be made to fit, especially given 
that it could be rearranged so that only one cpw-slotline transition was 
needed. In fact, by using radial stubs and at the cost of some defor- 
mation, one balun structure could be made to fit. However, as will 
be discussed later, because the radiation region of the spiral was wider 
than anticipated, 'several extra turns needed to be added to the spiral’s 
center. This greatly reduced the amount of space available at the cen- 
ter of the hollow spiral, making this implementation of an integrated 
planar feed impractical. 

• Balanced Twin- wire Feed 

The actual implementation of the balanced twin-wire feed is quite 
straightforward; some complexities exist, but techniques exist to over- 
come them [8]. On the other hand, the design and integration of the 
balun required to convert from the unbalanced coaxial line to the re- 
quired balanced twin-wire line is more complex. Several balun designs 
have been documented [4,8], but they either require a deeper cavity 
than allowed by the design specifications, or a relatively large external 
circuit. Hence, while the feed itself would work quite well, the idea was 
put aside until a more reasonable balun structure could be developed 
or acquired. 

• “Infinite-Balun” Feed 

Initially, because _of size constraints, there was no obvious way to im- 
plement this feed. The spiral arms are quite narrow, making it difficult 
to attach the coaxial lines to them without changing the driving-point 
impedance of the antenna significantly. To circumvent these problems, 
a slight modification was proposed, which amounted to replacing the 
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coaxial line with a microstrip line (this technique was later found in 
[10]). Because the width of a constant impedance microstrip line is 
determined by the substrate thickness, an “infinite-balun feed can be 
designed that has very little effect on the overall input impedance of 
the spiral antenna by choosing the appropriate substrate. 

Having decided to use the microstrip “infinite-balun” feed, there was 
some concern about the amount of coupling between the microstrip 
transmission line and the radiating slots (the microstrip line runs rela- 
tively close to the slots). However, by again choosing the appropriate 
substrate, an optimum design can be found which minimizes both the 
microstrip losses and the microstrip-slotline coupling. This is discussed 
further in a later section. 

For the microstrip implementation, there are two possible methods of 
exciting the slot at the antenna’s center. The first is similar to the 
method used with the coaxial “infinite-balun” feed — the microstrip is 
shorted to the opposite side of the slot, thus creating a current maxi- 
mum at the slot and maximizing the magnetic' field coupling between 
the microstrip and the slotline. This may be done with either a plated- 
through hole (via), or by cutting a slit, pulling the microstrip through, 
and soldering it to the other side, as shown in Figures 3a and 3b, re- 
spectively. 

The second method uses a A/4-long open-circuited stub, as shown in 
Figure 3c. The microstrip spirals inward towards the center of antenna 
as usual, crosses the slot, and terminates in an open-circuit approx- 
imately one-quarter of a center-frequency guide-wavelength away. It 
seems clear that each technique has its advantages and disadvantages. 
While the first will more accurately give a short circuit, and possibly 
better coupling and a better match, the via holes make it harder to 
fabricate. On the other hand, the second technique is considerably 
easier to build, but may not couple as well to the slot or give as good 
an impedance match. Moreover, it is relatively narrow band. For our 
design, since the bandwidth specification is quite narrow, the A/4-long 
open-circuit feed will be used. 


25 



6 Spiral Antenna Structure 

There are many radiating structures that are “frequency-independent.” In 
fact, any antenna that obeys the “angle principle,” as defined by Rumsy 
[7], will be frequency independent. However, to fit this criterion, most an- 
tenna structures also require infinite size. Clearly, a practical “frequency- 
independent” antenna must also obey the “truncation principle” [7], which 
states both that the current on the radiating structure must approach zero 
with distance away from the feed, and that, in addition, the radiation pattern 
of the structure approaches a limiting form with increasing frequency. Defin- 
ing the “active region” of the antenna as the area where the amplitude of the 
current diminishes rapidly due to radiation allows the “truncation principle” 
to be restated in a more useful way. It is seen to imply that the current 
on the antenna must greatly diminish in amplitude soon after it has passed 
the active region, so that the truncation of the antenna will not disturb it. 
Several well-known antennas obey these principles, notably the equiangu- 
lar spiral and the various forms of the log-periodic antenna. However, the 
archimedean spiral does not. It is interesting to note that it also doesn’t 
follow the “angle-principle,” and hence is not a true frequency-independent 
antenna, even though it appears so in practice. 

To truncate the archimedean spiral, even far beyond its active region, 
some sort of termination must be applied at the truncation point to minimize 
the reflections due to the abrupt change in impedance. To accomplish this, 
some idea of the size and position of the active region must be established 
so that the truncation is not applied too close, and the efficiency of the 
antenna reduced. An approximate method for determining the location and 
dimension of the active region is suggested in [8], and used to great effect in 
this study. Instead of plotting the current on the spiral arms, however, the 
field in the slot arms is plotted. Actually, a quantized version of the phase of 
the slot fields is plotted. In the plots shown in Figures 4 and 5a-c, the black 
lines indicate the phase <f) s.t. 0° < <f> < 180°, and the grey lines indicate 
the phase s.t. 180° < <f) < 360°. The active region can be easily identified 
by locating the areas where there are more that one or two lines of the same 
shade adjacent to one another. Figure 4 shows the fields in the spiral slots 
for the hollow spiral design. Comparison with Figure 5b, which shows an 
example of the standard spiral design, demonstrates the usefulness of this 
technique. It is readily seen that the active region for the spiral antenna at 
this frequency (the center frequency of the antenna) extends inward further 
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than do the slots in the hollow spiral design. Figures 5a-c show the same 
plots for the standard spiral design. Figure 5a shows the fields at the low 
frequency limit of 955 MHz, Figure 5b at the center frequency of 1100 MHz, 
and Figure 5c at the high frequency limit of 1247 MHz. These frequency 
limits were determined by keeping the bandwidth constant throughout the 
scaling procedure mentioned at the beginning of this report. It is interesting 
to verify that the active region moves away from the center as the frequency- 
decreases. 

From the the figures discussed above, it is apparent that the hollow' spi- 
ral design needs to be extended significantly inward towards the center, to 
accomodate the rest of the active region. However, this greatly reduces the 
hollow area in the center of the antenna, making it impossible to use the 
coax-cpw-slotline feed in its present form. For this reason, the hollow spiral 
design was set aside. 

Having determined the size and location of the active region, the termi- 
nation of the antenna can now be discussed. A good termination is critical 
to the operation of a spiral antenna because, if poorly executed, it is the 
greatest contributor to increased axial ratio. While -experiment has shown 
that the spiral arms can be truncated inside the active region [8], it has also 
shown that this requires very accurate terminations, and that a reduction 
in gain and pattern performance is to be expected. Although termination 
outside the active region is not as demanding, many of the same methods 
are used because they are simple and work well. Most of these involve the 
tapered application of a lossy material such as absorber foam or resistive 
paint to the spiral arms or slots. These can either be specially applied, or 
can be made of the foam in the absorbing cavity. There are many examples 
of these sorts of terminations in the literature [2,8,9]. For the antenna in 
question, a tapered application of lossy foam below the antenna will be used, 
as shown in Figure 6, because a cavity will not be used initially. 

Having qualified the fundamental issues involved in designing a spiral 
antenna, the actual dimensions can now be determined. The number of turns, 
the growth rate, the width of the slotline and of the microstrip feedline, and 
to a lesser extent, the overall diameter of the antenna, are all inter- related. 
The driving issue of the design eventually becomes the feed structure. The 
allowable amount of coupling between the microstrip line and the slotline and 
the thickness of the substrate effectively combine to determine the minimum 
allowable growth rate for the spiral. The thinner tlie substrate, the closer 
together the spiral slots can be for a given amount of coupling with the 
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microstrip feedhne, and the smaller the growth rate can be, yielding a higher 
performance antenna. Eventually this is limited by construction feasibility 
issues. The width of the slotline also affects the choice of the growth rate, but 
to a much lesser extent. The characteristic impedance of a wide slotline on 
a thin subtrate (slot width » substrate thickness) is virtually independent 
of the slotline width, allowing a wide range of choices. The number of turns 
is basically determined by the growth rate of the spiral, and the required 
operating bandwidth, as is the outer diameter of the antenna. 

For this antenna, the following parameters give a good compromise be- 
tween microstrip loss and coupling, as discussed above. For a dielectric thick- 
ness of 18 mils, the 50 0 microstrip feedline is 32 mils wide, and the 100 Q 
slotline is 44 mils wide. The growth rate was chosen to be a = .166, which 
spaces the slotlines 205 mils apart, center- to-center. Figures 5a-c were gener- 
ated using these values, and indicate that approximately 7 turns are required 
to adequately cover the required frequency band and allow for a good termi- 
nation. 


7 Experimental Results 

A spiral slot antenna was constructed according to the dimensions mentioned 
in the preceding section, and both the input impedance and pattern charac- 
teristics were measured. 

Figure 7 shows the reflection coefficient of the antenna, measured at the 
antenna’s connector, between 500 MHz and 5.0 GHz. For the frequency 
range of interest (950 to 1250 MHz), the reflection coefficient is approxi- 
mately — 17.5 dB, corresponding to a VSWR of 1.31. Contrary to the above 
discussion, the A/4-long tuning stub does not seem to restrict the bandwidth 
of the antenna significantly. An investigation of this effect is currently un- 
derway. 

The radiation patterns were measured for five different antenna configu- 
rations: 

1. The antenna, with no modifications. 

2. The antenna, with a ring of absorbing material placed on the slot side to 
terminate the slots, as discussed above. The ring was made of Eccosorb 
LS-26, 0.125" thick, with I.D. = 5.625" and O.D. = 9.00". The ring of 
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absorber was held against the antenna with a ring of thin cardboard, 
and fastened with eight nylon screws around its periphery. 

3. The antenna, with another absorber ring inside the first. The second 
ring was made of the same thickness material, with I.D. — 5.375", and 
O.D. = 5.625". 

4. The antenna, with the inner absorber ring removed, and a second ring, 
with the same dimensions as the first, placed on top of the first. Both 
rings were held on with the same structure used in configuration 2. 

5. The antenna, with both layers of absorber, and a ground plane placed 
on top of the absorber rings. All three layers were held on with the 
same structure as above. 

Patterns taken at the center of the frequency band of interest (1100 MHz) 
for each of the five modifications mentioned above are shown in Figures 8a-e, 
respectively. 

To determine the axial ratio of the radiation pattern, and thus the an- 
tenna’s ability to recieve or radiate circularly polarized fields, a rotating 
source antenna is usually used. This allows the polarization of the inci- 
dent field to be changed rapidly compared to the azimuthal rotation of the 
antenna-under-test. As the antenna-under-test rotates in azimuth, the vari- 
ation of the received level then indicates the dependence of the gain on the 
incoming field polarization at every azimuth angle, from which the axial ratio 
can be derived. 

Unfortunately, the facility where the radiation patterns in Figures 8a-e 
were measured is not equipped with a rotating source antenna, making this 
technique difficult at best. To give some idea of the axial ratio however, the 
following scheme was used. At boresight, the source antenna was first rotated 
to maximize the signal received by the antenna-under-test, and an azimuthal 
pattern taken. Then, again at boresight, the source antenna was oriented 
to minimize the received signal, and another azimuthal pattern taken. The 
main disadvantage of this technique is that it only gives the true axial ratio 
at boresight — the axial ratio at all other angles is only an optimistic estimate 
of the true axial ratio. For the patterns in Figures 8a-e, boresight is at 0°, 
and the two patterns are a result of the maximization and minimization of 
the received signal, as discussed above. 
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The radiation patterns for the unmodified antenna are shown in Figure 8a. 
At boresight, the axial ratio is approximately 6 dB, and at 180°, the axial 
ratio appears to be very close to zero. However, because the source polariza- 
tion wasn’t calibrated with the antenna-under-test at this orientation, this 
is an optimistic estimate at best. An investigation of the axial ratio around 
180° is currently underway. 

To reduce the axial ratio, as discussed above, a ring of absorber was used 
to create a tapered termination for the spiral slots. The radiation patterns 
resulting from this modification are shown in Figure 8b. The benefits of this 
modification, insofar as indicated by these patterns, are somewhat nebulous. 
The axial ratio at boresight did not improve, and the axial ratio at 180° 
appears to have deteriorated significantly. However, both the front-to-back 
and side-to-side symmetry of the patterns have greatly improved, and the 
location of the nulls is more realistic. In spite of the seeming decrease in 
pattern quality, these improvements to the symmetry give some credibility 
to the modification. 

In an effort to improve the axial ratio, another ring of absorber was added 
inside the first, as discussed above. If the slots weren’t adequately terminated 
previously, this should have improved matters, reducing the axial ratio some- 
what. However, as shown in Figure 8c, the axial ratio increased noticeably, 
and extreme pattern asymmetry resulted. While a better termination would 
not have these effects, an improperly terminated slot would reflect energy 
back along the slots, and cause similar problems. On examination of the 
second ring of absorber, a jagged edge was discovered (much rougher than 
that of the first ring). This would seem to bear out the necessity of a very 
precise termination, as discussed above. 

Still working under the assumption that the termination was not func- 
tioning adequately, the inner absorber ring was removed, and another ring 
of absorber placed on top of the first. The resulting pattern, shown in Fig- 
ure 8d, is virtually identical to the pattern in Figure 8b. This result would 
seem to indicate that, in fact, the slot termination was functioning properly, 
and that some other source must be responsible for the poor axial ratio. 
Further discussion of this follows in the next section. 

Accepting for the moment the poor axial ratio, a ground plane was placed 
behind the slot spiral to verify its unidirection operation. Figure 8e shows 
the resulting pattern. The axial ratio improved somewhat, and the rear 
radiation levels were reduced to an acceptable level. However, the overall gain 
of the antenna decreased by approximately 4 dB. Increasing the frequency, 
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however, verified that the ground plane had been placed too close to the 
slots, and had shorted some of the field in the slot. At 1500 and 2000 MHz, 
instead of decreasing the gain by 4 dB, the ground plane increased the gain 
by approximately 1.5 dB over the antenna without the ground plane. 


8 Future Work 

As discussed above, the poor axial ratio doesn’t seem to be a result of the 
reflections an improper slot termination would create. Based on the patterns 
shown in Figures 8a-e and on the fundamental geometrical symmetry of the 
antenna, another possible source of the high axial ratio is radiation from 
the microstrip feedline. The feedline’s radiation would have the opposite 
circular polarization as the spiral slots, greatly increasing the overall axial 
ratio. Feedline radiation would also cause a front-to-back asymmetry in the 
radiation patterns, as the microstrip would have a greater effect on the front 
of the antenna that on the back. As the ground plane behind the microstrip 
line is rather narrow, this mechanism seems rather likely. 

To eliminate this possibility, another antenna is currently being fabricated 
on a thinner substrate. This will not only place the microstrip line closer to 
the groundplane, but also make it narrower. Both of these modifications 
should serve to reduce the feedline radiation, and help to reduce the overall 
observed axial ratio of the antenna. 

Further investigation of the A/4-long tuning stub are also underway, with 
the dual goals of determining the actual workings of the feed and developing 
a narrow-band feed for application to this antenna. 
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Coax-cpw-slotline Baiun Feed - Fits 
into centec of Hollow Spiral Antenna 
(below) 






Figure 3a. Microstrip "Infinite Figure 3b. Microstrip "Infinite 
Baiun" Feed with Short- Baiun" Feed with Short- 

Circuit Via. Circuit Slit. 
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Figure 3c. Microstrip "Infinite Baiun" Feed 
with Xg/4 Open-Circuit Stub. 


35 



(cm) 

0° ^ <(> < 1 80° 180° < <j> < 360° 

Figure 4. Quantized Phase (<J>) in Hollow 
Archimedean Spiral Antenna 
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Figure 5a. Quantized Phase (<{)) at /= 955 MHz for 
Standard Archimedean Spiral Antenna 









Figure 6. Tapered Lossy Foam Absorber for 
Spiral Arm Truncation 
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Figure 7. Sll of Unmodified Spiral Antenna 
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Figure 8a. Azimuth Pattern of Unmodified 
Spiral Antenna 
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Figure 8c. Azimuth Pattern of Spiral with 
Inner Absorber Ring 
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Figure 8d. Azimuth Pattern of Spiral with 
Double-Thickness Absorber Ring 
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Abstract 

An efficient and accurate coax cable feed model is proposed for microstrip or cavity- 
backed patch antennas in the context of a hybrid finite element method (FEM). A TEM 
mode at the cavity-cable junction is assumed for the FEM truncation and system exci- 
tation. Of importance in this implementation is that the cavity unknowns are related 
to the model fields by enforcing an equi-potential condition rather than field continuity. 
This scheme proved quite accurate and may be applied to other decomposed systems 
as a connectivity constraint. Comparisons of our predictions with input impedance 
measurements are presented and demonstrate the substantially improved accuracy of 
the proposed model. 


I. Introduction 

The coax cable is widely used as a feeding structure for microstrip or cavity-backed patch 
antennas because of its simplicity and low spurious radiation. Indeed, abundant literature 
exists on the theoretical and experimental investigation of coax cable feeds [1] [3] . Most of 
these papers present integral equation techniques in conjunction with the pertinent Green’s 
functions. However, the Green’s function is only available for a certain class of geometries, 
and this limits the application of the integral techniques to those geometry designs. Also, the 
formulation must be modified and recoded for different antenna configurations corresponding 
to Green’s function variations. To avoid the complexity of the Green’s function, a hybrid 
finite element - boundary integral approach [5] was recently proposed to simulate cavity- 
backed antennas of arbitrary shapes. For antenna radiation it is observed that a simple 
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probe model with a constant current along the inner conductor linking the grounded base 
to the antenna element is straightforward and efficient. But the probe feed is only valid for 
thin substrates and this is consistent with the Moment Method (MM) results. Xo model 
an electrically thick substrate, in this communication we propose a more sophisticated feed 
modeling scheme in the context of the finite element method (FEM) using linear edge-based 
tetrahedral elements. In the next section, the formulation of the entire hybrid numerical 
system is first described for reference purposes. Xhe proposed feed model is then presented 
in section III on the basis of a TEM mode excitation. Model improvements are also discussed 
for the case when the TEM assumption at the cavity-cable junction does not hold. Input 
impedance comparisons between measurements and predictions based on this technique are 
presented in section IV. 


II. Hybrid FE-BI System 


The functional pertinent to the radiation by a cavity-backed antenna with a coax cable feed 
(as shown in figure 1) is given by 

F(E) = ^ JJJ V j(V x E) • ~“(V x E) - k^trE ■ E j dv 

-2kl JJ s ( E x z) • | JJ si ( E x z) ■ (l+ j^ VV ) G 0 (r,r')dS'j dS 

-jk 0 Zo JJ c (ExH)-zdS, (1) 

where V refers to the cavity volume and the surface S encompasses the cavity aperture 
excluding the portion occupied by the metallic antenna elements*, c r and /i r denote, respec- 
tively, the relative permittivity and permeability; ko is the free space wave number, Zo is 
the free space intrinsic wave impedance, I is the unit dyad, and Go(r,r') is the free space 
Green’s function with r and v f denoting the observation and integration points, the surface 
C is the cross section of the coax cable at the cavity-cable junction. 

Following the standard discretization procedure [5], we obtain a system of equations of 

the form 


N c N s e Nc € 

E {KH*?}} + E {isy{£*}} + E 


^ e dF c (E,H) 


e=l 


eC 5 




dE, 


= 0, 


(2) 


where the explicit expressions for Aij and Bij may be found in [5] and the functional term 
F c is the surface integral on C in (1) . It is apparent that the system excitation is associated 
with the functional Fc whose evaluation is the subject of this communication. 


III. Coax Cable Feed Modeling 

To proceed with the evaluation of 

F c = -jk 0 Z Q JJ(ExH)-z dS , (3) 
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a boundary constraint relating E to H is needed. To this end, we assume a TEM mode on 
C and consequently the fields within the cavity may be expressed as (see Figure 2) 


E = 2^' 1 + r > 7 ' 


(-0 


where t TC is the relative permittivity inside the coax cable; T denotes the reflection coefficient 
measured at 2 = 0 and I 0 is the given input current source at the same location. Also, (r, 4 > , 2 ) 
are the polar coordinates of a point in the cable with the center at r = 0. To simplify the 
analysis, we introduce the quantities 

— 5 ^ (1 + r) ' A “ = £ (1 - r) ' (5) 


Hence, 




( 6 ) 


and from (5) it follows 



( 7 ) 


which is the desired constraint at the cable junction in terms of the new quantities h 0 and 
e 0 - Note that eo and h 0 are field coefficients as new unknowns in place of the fields E and 
H, and it is therefore appropriate to rewrite Fc in terms of these new coefficients. To do so, 
we substitute (5) and (7) into (3) and upon making use of the axisymmetric field property 
we obtain 

F c = -2njk 0 Z 0 eoh s 0 rc ln(-), (8) 

a 


where a and b are the radii of the inner and outer cable conductors. The superscript src 
stands to indicate that h 0 is treated as a source term in the extremization of the functional. 

We choose the linear edge-based tetrahedral elements to discretize the cavity and the 
corresponding mesh on the cross section C is shown in Figure 2(b). In this formulation, the 
field across the p th edge, p=l,2,..., Nc {N c is the number of cavity mesh edges on C), is 
set to a constant as dictated by the linear edge-based expansion function inside the cavity. 
However, the cable TEM modal fields (4) behave as 1/r and this modeling inconsistency 
makes it difficult to apply the tangential field continuity condition at the cable junction 
( i.e. over the aperture C). To overcome these difficulties, we can relate the fields across 
the cable junction by recognizing that the potential difference between the inner and outer 
conductors must be the same as computed by the fields of the cavity or those in the cable 
region. Specifically, if the N* h edge of the cavity region borders and is also across the coax 
cable, from (4)— (6) it then follows that the appropriate equi-potential condition is 

A V = EAb — a) = e 0 ln~, i = N p (p = 1,2, ..., Nc)- (9) 

a 
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where Al' denotes the potential difference between the inner and outer surface of the cable. 
This condition simply provides a relation between the constant cavity edge field and the 
coax cable modal field. When used into the functional E c • it introduces the excitation into 
the hybrid finite element system without a need to extend the mesh inside the cable or to 
employ a fictitious current probe. It remains to rewrite F c in terms of i.e. the field value 
of the edges bordering the cable and to do so we substitute (i) and (9) into (b). Then upon 
taking into account all Nc cavity mesh edges on the cable junction, we obtain 


7T 


F c = --^jk 0 Z 0 {b- a) < 


fid y/^rc b 


7T 



( 10 ) 


In this expression, rather than representing the functional Fc in terms of a single edge field, 
we made use of the average field across the cable as computed by the totality of the equal 
element fields on the cable’s aperture (because of the axisymmetric property, all elements 
fields at the cable’s aperture are equal). The factor inside the curly brackets of (10), with the 
superscript src, functions as a source in the extremization process. Hence, the extremization 
of (10) yields 


dF c 

dE t 


--njk 0 Z 0 (b - a) 



x/tZcb-a „ \ 

Z 0 ln\ 7 



= UiEi-Vi, i = N p (p = 1,2, ...Nc), 

(11) 

where 

1 , H-ib-af 
U> - J-Kk 0y /e TC 

a 

(12) 


Vi = j^k 0 Z 0 (b- a)I 0 . 

(13) 


We observe that the ‘constant cavity field’ along each mesh edge at the cable junction is 
just a fictitious field representation and its meaningful physical interpretation is governed by 
the equi-potential constraint (9). To proceed, we assemble the FEM system together with 
(11). Specifically, each Ui is added to the N c diagonal entries of the finite element matrix 
which is associated with the Nc edges bordering the coax cable. Also, the excitation column 
of the hybrid system is nullified everywhere except for the Nc entries which are set to V t . 
Once the hybrid FE-BI system is solved [5], the input admittance at z = 0 is calculated from 


Y in = TT <f H • r rd<t> 

VO J 2ir 

2/p 1 

e 0 /n(£) Z c 


(14) 


where Z c is the characteristic impedance of the coax cable. 

In the above feed model we assumed the presence of only the dominant (TEM) mode at 
the cavity-cable junction, an assumption which may not be suitable for certain applications. 
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To overcome this limitation, one approach is to extend the mesh (say, a distance d) into 
the cable. The equi-potential condition will then be applied at z=-d, where all higher 
order modes vanish. This scheme requires a more suitable expansion for the fields in the 
— d < z < 0 section to avoid the complication of extending the tetrahedral mesh into the 
cable and, thus, retain the efficiency of the equi-potential feed model. Since the antenna is, 
in most cases, operated in the frequency range far below the cut-off of the first higher order 
mode of the coax cable, the field distribution near the junction C will still be dominated by 
the fundamental TEM mode. With this understanding, a possible suitable expansion for the 
field in the coax cable (using shell elements rather than tetrahedrals) is 






(15) 


where q=r, </> or z, i=l ,2,3 or 4 and N^(r, 4>, z) is defined as 


e. 


N* = 

Aq b A q, 


(<?6 - 9fr)(?c - q'cHa 


(16) 


with q a ,q b and q c representing r,d> and z in a cyclic rotation. Also, i denotes the edge 
number along each coordinate, and A q a is the width of the edge along the q a direction. The 
correspondence between the edge numbers and the node pairs for each coordinate^, cf> or z ) 
is given in Table 1 along with the definition of the primed parameters and e,- which carries 
the sign of the vector basis. 


TABLE1 
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q 
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Z\ + Az 

r 

T 

"Fay 

- 


4>x 

Z\ 4 Az 


3 

(8-7) 

+ 


4>i 

Z\ 


4 

(5-6) 

- 


4>\ 4- A 4> 

Z\ 


i 

(1-4) 

+ 

rx + Ar 


Z\ + Az 

<f> 

2 

w 

- 

ri + Ar 


Z\ 


Q 

(6-7) 

+ 

r\ 


Z 1 


El 

(2-3) 

- 

T\ 


z\ + Az 


i 

(1-5) 

+ 

n -f Ar 

4 > i + A <j> 



2 

(2-6) 

- 

ri 

(pi + A<f> 



3 

(3-7) 

+ 

n 

<h 



4 

(4-8) 

- 

r + Ar 

<h 




When an axisymmetric field property is assumed, the expansions reduce to 

N i _ ( g ~ fi) N 2 _ (f ~ ~ A *) 

T ~ Az ’ r “ Az 

M i _ (r-ri) 2 _ (r-n - Ar) 

Ar ’ Ar 


(17) 
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where (n.c!>i,2i) denote the coordinates of node 1. The accuracy of expansion (17) is demon- 
strated in Figure 3, where we show that only 3 elements are needed along the radial direction 
for the accurate prediction of the dominant field distribution. When compared to the con- 
ventional linear tetrahedral elements, the efficiency of this expansion is apparent (i.e., many 
more tetrahedrals are needed to model the same cable region). 


IV. Results and Conclusion 

To validate our proposed feed simulation, two circular patch antenna configurations were 
used for calculation. One patch antenna was of radius 1.3 cm printed atop of a circular 
cavity (radius=2.1 cm) filled with a dielectric (c r =2.9) material 0.41 cm deep. For this 
patch, the feed was placed 0.8 cm from the center and the input impedance was measured 
over the band 2-5 GHz. In Figure 4 we compare the measured input impedance with data 
computed on the basis of the proposed equi-potential feed model. Clearly, the results from 
measurements and the equi-potential model are in excellent agreement whereas the probe 
model yields substantially inaccurate results near resonance. 

Figure 5 shows the comparison between measurements and calculations for another patch 
antenna whose input impedance was measured by Aberle and Pozar [6]. This patch had a ra- 
dius of 2.0 cm and the 0.218 cm thick substrate had e r =2.33 and a loss tangent tan6=0.0012. 
The feed was located 0.7 cm from the center, and for our FE-BI calculation the patch was 
placed in a circular cavity of 2.44 cm in radius. As shown in Figure 5, the equi-potential 
model is again in excellent agreement with measurements, as opposed to the results by the 
probe model in [6] . 

In conclusion, the presented equi-potential feed model has been shown to be extremely 
accurate in modeling coax feed structures. Moreover, its implementation in the context of a 
finite element formulation is very simple and as easy to implement as the probe feed. It was 
also demonstrated how the proposed feed model can be generalized to the case of asymmetric 
feed structures where evanescent modes may be present. 
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having the following specifications: patch radius r=13mm; cavity radius R=21.1mm; 
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part; (b) Imaginary part. 
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Figure 1: Illustration of a cavity-backed patch antenna with a coax cable feed. 



(a) 


(b) 


Figure 2: (a) Side view of a cavity-backed antenna with a coax cable feed; (b) Illustration 
of the FEM mesh at the cavity-cable junction. 


55 


real part of input impedance (Ohms) 


Comparison of E-FiekJ Along (1 -Wavelength) Coax Cable 



segment number 
(a) 


Comparison of E- Field Distribution Along R 



(b) 


Figure 3: Field distribution in a shorted coax cable as computed by the finite element method 
using the expansion (14). — : analytical; xxx: numerical, (a) Field coefficient e 0 along the 
length of the cable (leftmost point is the location of the short); (b) Field along the radial 
coordinate calculated at a distance A/4 from the shorted termination. 



frequency (GHz) 

(a) 



Figure 4: Measured and calculated input impedance for a cavity-backed circular patch an- 
tenna having the following specifications: patch radius r=13mm; cavity radius R=21.1mm; 
substrate thickness t=4.1mm; e r =2.4; and feed location :r/=0.8 cm distance from center. 
Results based on the simple probe model are also shown for comparison, (a) Real part; (b) 
Imaginary part. 
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(a) (b) 


Figure 5: Measured and calculated input impedance for a circular patch antenna [6] having 
the following specifications: patch radius r=2cm; substrate thickness d 0.21844cm, feed 
location from center *,=0.7011; e r =2.33; i<m6=0.0012. [6]. — : measurement; xxx: this 
method; o o o: probe model [6] (a) Real part; (b) Imaginary part. 
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Abstract 

A dual-stacked patch antenna is analyzed using a hybrid finite 
element — boundary integral (FE— BI) method. The metallic patches 
of the antenna are modeled as perfectly electric conducting (PEC) 
plates stacked on top of two different dielectric layers. The antenna 
patches may be of any shape and the lower patch is fed by a coaxial 
cable from underneath the ground plane or by an aperture coupled 
microstrip line. The ability of the hybrid FEM technique for the 
stacked patch antenna characterization will be stressed, and the EM 
coupling mechanism is also discussed with the aid of the computed 
near field patterns around the patches. 


I. Introduction 

Since being proposed, the dual stacked patch antennas received tremendous 
attention mainly because of the drastic improvement of their operational 
bandwidth in comparison to the inherent narrow band feature of single patch 
configurations. Further bandwidth improvement was also reported by em- 
ploying an aperture coupled microstrip line, where the lower patch is electro- 
magnetically coupled to a microstrip line through an opening at the ground 
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plane. This design was shown capable of achieving over 10% bandwidth. 
Certain numerical approaches for extensive analysis of non-rectangular sin- 
gle patch antennas with an aperture coupled feed have been investigated 
(e.g. [1, 2]). Unfortunately reports for non-rectangular stacked structure 
analyses have not been seen in the literature, except for experimental results 
published in [3, 4]). An analysis was only reported in [5] where the Hankel 
transform was employed on the assumption that the feed structure is absent 
from the analytical system. This implies that any change in the resonance 
due to the feed can not be modeled using this technique. A comprehensive 
simulation for stacked non-rectangular patch configurations is therefore de- 
sired for design purposes. This is the motivation of this report and we will 
discuss both probe and aperture coupled feeds. The former offers flexibility 
in the feeding network layout and makes the feed design separable from that 
of antenna elements [3], whereas the latter design, while challenging for most 
of the existing methodologies, is readily incorporated into our proposed hy- 
brid FEM formulation, especially when the patches and the coupling opening 
(aperture) are non-rectangular for wider bandwidth designs. 

In the following sections, we first briefly describe the hybrid FEM formu- 
lation as applied to the patch antenna analysis. More detailed discussions 
of the numerical approach may be found in [6] and [7]. Using our numerical 
simulator for characterizing the antenna’s dielectric and geometric parame- 
ters, we concentrate on the microstrip feed line modeling with the aid of an 
FEM connectivity scheme. Some results for a circular stacked-patch antenna 
will be presented to show the ability of the technique, and typical computed 
3-D near field distributions around the patches inside the two dielectric layers 
are displayed for visualization which reveals the insight of the EM coupling 
mechanism between the two patches. 


II. Hybrid FEM Formulation 

Consider the configuration of a stacked circular patch antenna shown in Fig- 
ure 1, where a cavity encompassing the patches is recessed in the ground plane 
to simulate the truncation of the computational region V. As described in 
[6], the finite element method (FEM) is employed to model the fields in the 
three dimensional volume, whereas a boundary integral equation on the aper- 
ture of the cavity (excluding the metallic patch area) is built up based upon 
the equivalence principle to associate the radiated fields in the region above 
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the ground plane to those in the cavity so that the continuity conditions of 
the tangential fields across the aperture is inherently imposed. (Interested 
Readers are referred to [6] for detailed discussions and mathematical deriva- 
tions.) It should be noted that the edge-based tetrahedral elements were 
used for the 3-D cavity fields and the triangular surface elements for the 
aperture fields simulation. The reasons are twofold (a) the edge elements 
have shown a precise representation of the vectorial electromagnetic fields 
without spurious modes; (b). arbitrarily shaped geometries and more so- 
phisticated feeding schemes (e.g. aperture coupled stripline, microstrip line, 
equi-potential model for coaxial cable [8], etc.) are easier to model using 
edge elements. After the standard discretization and assembly procedures, 
we obtain a system of equations in a matrix form as 

Ny Ns 

^ [Aivj] {£j} + ^ = {A,-} , * = (1) 

j i 

where i v and i s represent the unknown indices in the cavity and on the 
aperture, respectively, and the former contains the latter in global numbering 
since the aperture fields E * are part of the cavity fields Ej . Nv and Ns are 
the number of cavity and aperture edges that do not lie on any PEC surfaces. 
It is apparent that the total number of the system unknowns is N v and, in 
most cases, Ns < N v . This property is important because [A,*y] is a sparse 
system although [B^j] is a full matrix. Thus the storage requirement is 
insignificant for this exact boundary truncation, which governs the modeling 
accuracy for the radiation problems. 

In the case of stacked patch antennas, the lower patch is usually placed 
close to the cavity’s base, which is treated as a grounded PEC wall. Hence, 
inside the FEM cavity, the lower patch should also be modeled as a PEC 
plate. Two separate dielectric layers above and below this patch must be 
considered to account for the substrate effects. Indeed, an attractive fea- 
ture of the FEM over integral equation method is its capability to model 
multilayer structures with the same , case as it does with the homogeneous 
structures. 
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III. Feed Simulations 


Probe Feed 

For the radiation analysis of a probe-fed patch antenna, we use the infinitesi- 
mally thin current source which is chosen to coincide with some mesh element 
edge. In this fashion, the excitation term {At"} is usually modeled such that 
all the entries are set to zero except for those coincident with the probe(s) 
which link the patch and the ground plane. Specifically, if the probes are 
assumed to have a constant current 7„i = 1,2, M v (M p is the number of the 
probes), then the excitation entries associated with the corresponding edges 
may be written as 

Ki = - jkZoUU , ( 2 ) 

where k is the free space wave number, Z 0 is the free space intrinsic wave 
impedance and /, is the length of the i th edge. The phase difference between 
different probes may be readily represented by the phasor notation of each 
complex current /,-. Once the system of equations are solved, the input 
impedance is found to be 


where Ei denotes the calculated field along the i th edge of length The 
radiation pattern may be obtained by integrating the computed aperture 
fields with the far-zone approximation of the dyadic Green’s function in the 
presence of an infinitely large PEC ground plane. 


Microstrip line Feed 

When the stacked patch antenna is fed with a microstrip line network under- 
neath the ground plane (cavity’s base) via a coupling aperture (see Figure 2), 
special treatment for the feed structure must be considered in the FEM for- 
mulation. This is because the microstrip line is usually designed to have 
different size and shape as compared to the cavity’s geometric parameters. 
In addition, the microstrip line used as a transmission line, in most cases, 
carries traveling waves, thus, changing the field behavior in the cavity. Hence, 
the conventional FEM simulation treating the entire 3-D domain in a same 
fashion is not appropriate for this problem. 
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Consider the structure to be modeled as shown in Figure 3. In this case, 
it is appropriate to separate the computational domains because of the small 
element size required in modeling the feed structure. One difficulty encoun- 
tered when this decomposition is implemented is how one must model the 
coupling through the aperture. (We consider in this report a rectangular 
aperture which has been extensively employed in practice. Other aperture 
shapes for wideband consideration have also been seen in the literature and 
this formulation can be readily extended for those designs.) As described in 
[6], the cavity fields are discretized using tetrahedral elements, whereas in 
the microstrip line rectangular brick elements are considered the best candi- 
dates for field modeling. The selection of brick elements is particularly useful 
because of the geometry and the field property of the microstrip line, and it 
is shown that this choice retains better accuracy and efficiency in comparison 
to other element selections. As shown in Figure 4, although both types of 
elements employ edge-based field expansions, the meshes across the common 
area (coupling aperture) are different, and this causes the inconsistency for 
the field representation in the formulation. In the implementation, this im- 
plies that the usual boundary condition for the tangential field continuity 
across the aperture can not be enforced in a straightforward manner. 

We observe that the aperture is narrow in shape along the microstrip line 
and this fact allows us to assume a ‘static’ field feature for a given frequency. 
Therefore, the ‘potential’ concept may be applied to the problem in order to 
relate the fields below and above the aperture. Specifically, we enforce the 
‘equi-potential’ continuity condition, which states that: the potential across 
the narrow aperture due to the field distribution in the cavity region must be 
equal to that due to the microstrip line excitation. If the edges are classified 
as 


in mesh (b): 

• Ef j = 1,2,3, 

• Ef j = 1,2,3, 
in mesh (c): 

• Ej j = 1,2,3, .. 


vertical edges 
diagonal edges 

vertical edgesonly 


then the ‘equi-potential’ continuity condition is expressed as 


Ef = 




( 4 ) 
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(5) 


Ef 


-( Cj E' + t j+1 E' +1 ), 


where 

«, = («) 

t and d are the lengths of vertical and diagonal edges, respectively. (Thus, t 
happens to be the width of the narrow rectangular aperture.) 

The connectivity scheme for entirely different computational domains 
may be extended by generalizing this concept. It is apparent that this ap- 
proach makes the FEM implementation straightforward for different geom- 
etry/size domains that would be significantly inefficient if only one type of 
elements were used for modeling the structure. In addition, this technique 
ensures the system condition since the number of distorted elements in the 
mesh are minimized. 


Results 

As an example, a circular stacked patch antenna was used as shown in Fig- 
ure 1, where the geometry parameters were designed for the antenna to be 
operated in the frequency range from 7 to 9 GHz. When excited using a 
microstrip line, it is observed that the field across the coupling aperture has 
a Gaussian-like distribution, as expected. Figure 5 shows the gain com- 
puted as a function of frequency and as seen, the design is indeed wideband, 
as promised. Figure 6—9 show the radiation patterns for both polarizations 
in the E-plane and H-plane, respectively, and Figure 10 plots the input 
impedance, whose VSWR is 2:1, corresponding to a bandwidth of 15%. In 
Figure 1 1 , we show the magnitude of the computed field distribution on the 
antenna’s radiating aperture when a probe feed was employed for excitation. 
It is noted that the energy is concentrated on the two opposite sides of the 
top patch and that the PEC cavity wall has little effect on the field near 
the antenna edges. Also, we observe from Figure 12 that the field couples to 
the lower patch from two locations even though only one probe is used for 
excitation. 
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Figure 1: Illustration of A Cavity-backed Stacked Circular Patch Antenna, 
which could be fed by a probe or by an aperture coupled microstrip line. 
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Tj = 0.65 cm 
r 2 = 0.75cm 
t = 0.08128 cm 
d-t = 0.18128 cm 
d= 0.26256 cm 
s= 0.2 145cm 

D = 2.6 cm* * Calculations were done using this cavity size. 

Thus, performance will remain unchanged 
when array elements are separated by comparable 
distances. Tighter placement requires further 
investigation to maintain performance within specs. 

Figure 2: A stacked circular patch antenna with a microstrip line feed. The 
design of geometric and dielectric parameters will be used for modeling. 


slotwidth=0.04 cm 
Xg- guided wavelength 
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Antenna Elments 



Figure 3: Cross-section of an aperture coupled patch antenna, showing the 
cavity region I and the microstrip line region II for two different FEM com- 
putation domains. 
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(a) 



Figure 4: (a), a rectangular aperture; (b). a typical mesh from cavity region 
I ; (c). the uniform mesh from microstrip line region II. 





Gain VS Frequency for a Stacked Circular Patch Anttena 



Figure 5: The gain pattern vs. freqency in GHz for circular stacked patch 
antenna. 
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Figure 11: The field magnitude distribution on the antenna’s radiating aper- 
ture. It is shown that the energy is concentrated on two opposite sides of the 
top patch and that the PEC cavity wall has little effect on the field. 
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Figure 12: The field magnitude distribution within the cavity of the antenna 
excited by a probe. It is shown that the energy couples to the lower patch 
from two locations even though only one probe is utilized for excitation. 



N95- 19801 


3 73 * 

P- > 

Performance of an Artificial Absorber 
for Truncating FEM Meshes 

Jian Gong and John L. Volakis 
Radiation Laboratory 

Department of Electrical Engin. and Comp. Sci. 

University of Michigan 
Ann Arbor, MI 48109-2212 


Abstract 

We investigate the effectiveness of an artificial absorber for truncating 
finite element (FE) meshes. Specifically, we present the implementation of 
a novel mesh truncation approach using a perfectly matched anisotropic 
absorber for waveguides and stripline circuits. This truncation scheme is 
useful in many applications, including antennas, scattering and microwave 
circuits. 


I. Introduction 

One of the most important aspects of finite difference and finite element imple- 
mentations is the truncation of the computational volume (see Figure 1). An ideal 
truncation scheme must ensure that the outgoing wave is not reflected backwards 
at the mesh termination surface, i.e. the mesh truncation scheme must simulate 
a surface which actually does not exist. To date, a variety of non— reflecting or 
absorbing boundary conditions (ABCs) have been employed for truncating the 
computational volume at some distance from the radiating or scattering surface. 
These conditions have also been extensively used in modeling microwave circuits 
and devices. The ABCs are typically second or higher order boundary conditions 
and are applied at the mesh termination surface to truncate the. computational 
volume as required by any PDE solution. Many 2D and 3D variations of ABCs are 
nowadays available and several successful applications of these have been reported 
in the literature. Among them, a class of ABCs is based on the one-way wave 
equation/method [1, 2, 3] and another is derived starting with the Wilcox Expan- 
sion [4, 5]. Also, higher order ABCs using the Higdon [6] or Liao [7] formulation 
and problem specific numerical ABCs have been successfully used particularly for 
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truncating needs in guided structures. There are several difficulties with tradi- 
tional ABCs. Among them is the need to apply them at some distance from the 
perturbance (leading to a large computational volume) and their limited accuracy 
particularly for truncating meshes atop of radiating elements or in guided struc- 
tures. The latter prompted the use of higher order ABCs using the Higdon or Liao 
formulation (see [8, 9, 10]) and a variety of numerical ABCs [11, 12]. Nevertheless, 
regardless of the improvements realized by higher order and numerical ABCs, they 
substantially complicate the implementation of PDE techniques especially when 
inhomogeneity is present and when parallelization of the codes is needed. 

An alternative to the traditional ABCs is to employ an artificial absorber for 
mesh truncation. Basically, instead of an ABC, a thin layer of absorbing material 
with c T = n r is used to truncate the mesh as illustrated in Figure 2. Typically, 
the absorber must be 0.15 A thick and so far only homogeneous layers have been 
investigated [13]. Such absorbing layers are therefore necessarily narrow band 
in addition to exhibiting a non-vanishing reflection coefficients for angles near 
grazing. One approach to overcome both of these difficulties is to employ a uniaxial 
rather than an isotropic material layer and in this communication we examine the 
performance of such a mesh truncation scheme in the context of a finite element 
implementation. More specifically, the uniaxial absorbing layer is employed for 
truncating meshes in waveguide and stripline propagation problems. 


II. Formulation 

Consider the rectangular waveguide and a shielded stripline structure in Figure 3. 
We are interested in modeling the wave propagation in these structures using the 
finite element method. For a general anisotropic medium, the functional to be 
minimized is 

7= f V x E • (p; 1 ■ V x E) - k% • E • EdV 
j v 

-f E x (ji~ l ■ V x E) • dS, (1) 

J *Sj n +Sout 

in which ju r and Z T denote the permeability and permittivity tensors whereas E 
is the total electric field in the medium. The surface integrals over Si n and S out 
must be evaluated by introducing an independent boundary condition and the 
ABC serves for this purpose but alternatively an absorbing layer may be used as 
illustrated in Figure 2. One approach in evaluating the performance of an absorbing 
layer for terminating the FE mesh is to extract the reflection coefficient computed 
in the presence of the absorbing layer used to terminate the guided structure. In 
this study we consider the performance of a thin uniaxial layer for terminating the 
FE mesh in a rectangular waveguide and a stripline. Such a uniaxial layer was 
proposed by Sacks et.al. [14] who considered the plane wave reflection from an 
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anisotropic interface (see Figure 4 ). If [i r and e r are tensors of the form 



fa 0 

0 \ 

Vr = = 

0 b 

° 


\ 0 0 

He) 

the TE and TM reflection coefficients at the interface 

R TE = 

cosOi — 

\feeos6t 

I — 


cosQi + 

V \cos9t 

R™ = 

y/icasB 

i — cosOi 

rr 


cosOi + 

\fccos9 t 


( 2 ) 


( 3 ) 


and by choosing a = b and c = £ it follows that R TE = R™ = 0 for all incidence 
angles, implying a perfectly matched material interface. Moreover, by setting 
a — a — j(3 = \— jf3 the reflected field of such a metal-backed uniaxial layer is 

|#(0,)| = e - 2 ^ otcos5 ' (4) 


where t is the thickness of the layer and 9{ is the plane wave incidence angle. Ba- 
sically, the proposed metal backed uniaxial layer can have a reflectivity of as much 
as -17 dB if / StcosOi = 0.16A o , where A 0 is the free space wavelength. The reflection 
coefficient (4) can be reduced further by backing the layer with an ABC rather 
than a PEC. However, the PEC backing is more attractive because it eliminates 
altogether the integrals over the surfaces. 

Below we present a number of results which show the performance of the pro- 
posed uniaxial absorbing layer as a function of the parameter (3 , the layer thickness 
t and the frequency for a waveguide and a stripline under single mode operations. 


Results 

Let us first consider the rectangular waveguide shown in Figure 3. The guide’s 
cross-section is chosen to propagate only the T E w mode and is excited by an elec- 
tric probe at the left. Figure 5 shows the mode field strength inside the waveguide 
which has been terminated by a perfectly matched uniaxial layer having (3 = 4.2 
and thickness t = 1.0 cm. Clearly, the field amplitude is fairly constant once away 
from the probe indicating a small VSWR. Plots which illustrate the performance 
of this uniaxial layer as a function of frequency for different values of (3 are shown 
in Figure 6. The plots show that for (3 = 4.2 the reflection coefficient is less than 
-17 dB over a broad range of frequencies. Interestingly, different values of (3 can 
have a better performance (see Figure 9) but most importantly, in all cases, the 
broadband performance is maintained. 

The performance of the perfectly matched uniaxial layer in absorbing the 
shielded stripline mode is illustrated in Figure 7. In this case, the stripline is 
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extended into the absorbing layer up to the metallic wall. We observe that the 
predicted reflection coefficient is now around -13 dB and nearly independent of h 3. 
In comparison to the rectangular waveguide example, there is a 4 dB reduction 
in absorption and this is likely due to the extension of the substrate within the 
absorbing layer. 

A most surprising observation on the performance of this uniaxial layer for 
truncating finite element meshes is extracted from Figure 8. In this figure, we plot 
the reflection coefficient as a function of layer thickness (actually as a function of 
the number of elements used in modeling the layer). It is seen that the performance 
of the layer for truncating the FE mesh is maintained reasonably well even when 
two to three sample elements are used within the layer. That is, the proposed 
perfectly matched uniaxial layer bares no computational burden on the simulation. 
In fact, it is even simpler to apply this anisotropic absorbing layer than the first 
order ABC since the PEC backing eliminates the surface integral from the FE 
formulation. 

Conclusions 

We examined the performance of a new perfectly matched absorbing layer for trun- 
cating FE meshes. Examples relating propagation in a rectangular waveguide 
and a stripline were considered and it was demonstrated that a very thin metal- 
backed anisotropic layer was sufficient to achieve a reflection coefficient which in 
most cases was -15 dB or less. Surprisingly, the layer could be modeled with as 
few as two elements across its thickness and its overall implementation was rather 
simple because of the metal backing. Such an absorbing layer is particularly useful 
in truncating meshes in guided structures and should be found as useful in the 
application of finite element methods to scattering prediction and open microwave 
circuit simulation. 
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Mesh Truncation Surface 



Figure 1: Mesh truncation using an artificial absorbing surface. 



Figure 2: Mesh truncation using a thin anisotropic layer. 
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E lectric 



(a), waveguide 


8x = £y= e z = 1-jp 

Mx = My = Ml- HP 

d + t=40cm t=1 cm 



d + 1 « 12.0 cm 
L = 2.38 cm 
h = 0.21 cm 


H * 1.06 cm 
w * 0.47 cm 


(b). Stripline 

Figure 3: A rectangular waveguide (a) and a stnpline (b) truncated using the 
perfectly matched uniaxial absorbing layer. 



Figure 4: Illustration of plane wave incidence on an anisotropic interface. 
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Electric 



E * = e y= e z= "HP 
m = n y = m*- i-jP 


d + t=40cm t= 1 cm 


Figure 5: Typical field values of the TE\o mode inside a rectangular waveguide 
terminated by a perfectly matched uniaxial layer. 



Figure 6: Reflection coefficient vs frequency for the perfectly matched uniaxial 
layer used to terminate the waveguide shown in Figure 3. 
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Figure 7: Reflection coefficient vs frequency for the shielded stripline terminated 
by the perfectly matched uniaxial layer. 



Figure 8: Reflection coefficient vs thickness t of the perfectly matched uniaxial 
layer for f = 4.2 GHz and j3= 4.2. 





Figure 9: Reflection coefficient vs (3 of the perfectly matched uniaxial layer for f 
= 4.2 GHz and t = 10 samples. 
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